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Abstract: We clarify some peculiar aspects of the perturbative expansion around a classi- 
cal fuzzy-sphere solution in matrix models with a cubic term. While the effective action in 
the large- N limit is saturated at the one-loop level, we find that the "one-loop dominance" 
does not hold for generic observables due to one-particle reducible diagrams. However, we 
may exploit the one-loop dominance for the effective action and obtain various observables 
to all orders from one-loop calculation by simply shifting the center of expansion to the 
"quantum solution" , which extremizes the effective action. We confirm the validity of this 
method by comparison with the direct two-loop calculation and with Monte Carlo results 
in the 3d Yang-Mills-Chern-Simons matrix model. From the all order result we find that 
the perturbative expansion has a finite radius of convergence. 
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1. Introduction 

Fuzzy spheres [1], which are simple compact noncommutative manifolds, have been stud- 
ied extensively. There are various motivations for studying the fuzzy spheres. First it is 
expected that the noncommutative geometry provides a crucial link to string theory and 
quantum gravity. Indeed the Yang-Mills theory on noncommutative geometry is shown to 
emerge from a certain low-energy limit of string theory [2]. There is also an independent 
observation that the space-time uncertainty relation, which is naturally realized by non- 
commutative geometry, can be derived from some general assumptions on the underlying 
theory of quantum gravity [3] . One may also use fuzzy spheres as a regularization scheme 
alternative to the lattice regularization [4]. Unlike the lattice, fuzzy spheres preserve the 
continuous symmetries of the space-time considered, and the well-known problem of chiral 
symmetry [5-19] and supersymmetry in lattice theories may become easier to overcome. 

As expected from the so-called Myers effect [20] in string theory, fuzzy spheres appear 
as classical solutions in matrix models with a Chern-Simons term [21-24]. The properties of 
the fuzzy spheres in matrix models have been studied in refs. [25-32]. One can actually use 
matrix models to define a regularized field theory on fuzzy spheres. Such an approach has 
been successful in the case of noncommutative torus [33], where nonperturbative studies 
have produced various important results [12,16,34]. These matrix models belong to the 
class of so-called large-iV reduced models, which are believed to provide a constructive 
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definition of superstring and M theories [35-37]. For instance the IIB matrix model [36], 
which can be obtained by dimensional reduction of ten-dimensional M = 1 super- Yang- 
Mills theory, is proposed as a constructive definition of type IIB superstring theory. In 
this model the space-time is represented by the eigenvalues of bosonic matrices, and hence 
treated as a dynamical object. The dynamical generation of four-dimensional space-time 
has been discussed in refs. [38-51]. 

In ref. [52] we have performed a first nonperturbative study on the dynamical properties 
of the fuzzy spheres, which appear in a simple matrix model. The model can be obtained by 
dimensional reduction of 3d Yang-Mills-Chern-Simons (YMCS) theory, and it incorporates 
various fuzzy S 2 solutions. The most important non-perturbative result was that there 
exists a first-order phase transition as we vary the coefficient of the Chern-Simons term 
(a). In the small-a phase, the effect of the Chern-Simons term is negligible and the model 
behaves almost as the pure Yang-Mills model (a = 0) [53]. In the large-a phase, on the 
other hand, a single fuzzy sphere appears dynamically 1 . In this phase Monte Carlo data 
agree very well with the one-loop results, which led us to speculate that the one-loop 
contribution dominates the quantum correction in the large- N limit [52]. 

In the present paper we first address this issue by direct two-loop calculation. While 
the "one-loop dominance" does hold for the effective action [27,31,32], we find that this is 
not the case for generic observables. The higher-loop contribution that survives the large- 
N limit actually comes from one-particle reducible diagrams, which do not appear in the 
calculation of the effective action. However, we may exploit the one-loop dominance for 
the effective action and obtain various observables to all orders from one-loop calculation 
by simply shifting the center of expansion to the "quantum solution" , which extremizes the 
effective action. We confirm the validity of this method proposed by Kitazawa et al. [31] 
by comparison with the direct two-loop calculation and with Monte Carlo results. 

From the all order result we find that the perturbative expansion has a finite radius of 
convergence, and the lower critical point of the first-order phase transition lies precisely on 
the convergence circle. We also reconsider the issue of the dynamical gauge group, which 
was previously discussed at the one-loop level [52]. The all order calculation of the free 
energy for k coinciding fuzzy spheres allows us to obtain a more definite conclusion on this 
issue. 

This paper is organized as follows. In section 2 we define the model and briefly re- 
view the results obtained in our previous work. In section 3 we perform explicit two-loop 
calculation of an observable and demonstrate that the two-loop contribution survives the 
large- N limit. In section 4 we obtain an all order result from one-loop calculation by 
shifting the center of expansion to the "quantum solution" . In section 5 we extend the all 
order calculation to more general observables. In section 6 we compare the results of the 
all order calculation with our Monte Carlo results. In section 7 we address the issue of the 
dynamical gauge group using the all order result for the free energy. Section 8 is devoted 

1 This work has been extended to matrix models which incorporate four- dimensional fuzzy manifolds as 
classical solutions [54-56]. While the fuzzy S 4 turned out to be always unstable, the fuzzy CP 2 and the 
fuzzy S 2 x S 2 can be stabilized in the large- N limit. 
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to a summary and concluding remarks. In appendix A we give the details of the two-loop 
calculation. 



2. Brief review of the model 

The model we study in this paper is given by the action 

S[A] = Ntr ^[^,4f + ^ae F/) V^j , (2.1) 

where Ap (fi = 1,2,3) are N x N traceless hermitian matrices. Here and henceforth we 
assume that summation is taken over repeated Greek indices. The classical equation of 
motion is given by 

[Ap,[Ap,A u ]]+iae ua p[A a ,Ap] =0 . (2.2) 

As a solution, we consider 

^=l/=al,, (2.3) 
where Lp is the ./V-dimensional irreducible representation of the SU(2) Lie algebra 

[Lp, L„] = it^vpLp . (2.4) 
This solution corresponds to the fuzzy S 2 , and since 

(Xpf = \ c? (N 2 - 1) 1 N , (2.5) 

the radius of the sphere is given by R = \a^/N 2 - 1. 

In ref. [52] we found that the system undergoes a first-order phase transition as we 
vary a. In the large-a regime the dominant configurations are close to the fuzzy sphere 
(2.3), while in the small-a regime the large- N property is similar to the pure Yang-Mills 
model (a = 0) [53], and the geometry of the dominant configurations is given by that of a 
solid ball. The two phases are called the "fuzzy sphere phase" and the "Yang-Mills phase" , 
respectively. When we discuss the large-iV limit in the fuzzy sphere phase, the natural 
parameter to fix turned out to be 

a = aVN. (2.6) 
The lower critical point obtained from Monte Carlo simulation is 

«cr ^ 2.1 , (2.7) 
which was reproduced later from the one-loop effective action as 2 



/ 512 

a cr = {J— = 2.0867794- •• . (2.8) 

In the fuzzy sphere phase, various observables agree well with the one-loop calculation. We 
therefore speculated that the one-loop dominance, which was previously claimed for the 
effective action [27], holds also for ob servables. 

2 This analytical result was informed to us by D. O'Connor after J.N. gave a seminar on the Monte Carlo 
results including (2.7) at the Dublin Institute for Advanced Studies (DIAS). Its derivation in section 4 is 
due to Y. Kitazawa (private communication). 
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3. Explicit two-loop calculation 

In order to see whether the one-loop dominance holds also for observables, we perform 

explicit two-loop calculation around the fuzzy sphere solution. We decompose into the 
classical background = aL^ and the fluctuation as 

A ll = X li + A li . (3.1) 

We introduce the following gauge fixing term and the corresponding ghost term 

S s .i. = -\Ntv([X„A,} 2 ) , (3.2) 

S f ^ = -Ntr([X li ,c][A li ,c]) . (3.3) 

The total action can be written as 

<Stotai = S + 5 g .f. + S gh (3.4) 

= S[X] + S kin + S int , (3.5) 

where the kinetic term 5km an d the interaction term S m t are given by 

S k in = \ N ^ (M X x, [X X , \]}) + Nti (c[X x , [X x , c]]) , (3.6) 

+-iaNe lwp tT(A li A v A p )-Nix([X li ,c][A li ,c]) . (3.7) 
The free energy W defined by 

e~ w = J dAdcdce- St °^ (3.8) 

can be calculated as a perturbative expansion 

oo 

W = J2Wj, (3.9) 

3=0 

where Wj represents the j-th order contribution. The first two terms are obtained as [21,52] 

a 4 

W = S[X] = --(N 2 -l) , (3.10) 



1 r i 

Wl = 2 ^( 2/ + 1 ) 1 °S a 2 1(1 + 1) . (3.11 



l=i 

In order to calculate W2, we have to evaluate the two-loop diagrams 3 depicted in figure 
1. The solid line and the dashed line represent the propagators of A^ and the ghost, 
respectively. The three-point vertices with (without) a dot represent the third (second) 
term in (3.7). 

3 The diagrams (a)~(d) are the same as the ones that appear in ref. [27]. The diagram (e) of ref. [27], 
which involves a fermion loop, does not appear in the present bosonic model. 
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Figure 1: The two- loop diagrams that appear in the perturbative calculation of the free energy. 
The diagrams (a) ~ (d) are one-particle irreducible, while the rest are not. 



As a fundamental observable, let us consider the vacuum expectation value of the 
action (S). The perturbative expansion of this quantity can be readily obtained from the 
results for the free energy (3.9). Let us define the rescaled action 



S(X, a) = A iVtr (-^[A^ A u f + a e^ p A^A u A^j 



and the corresponding free energy 

e~ w ^ = j dAdcdce- s ^ , 
which is related to the original free energy W = W(l,a) through 

W(X,a) = -(N 2 -1) logA + W(l,A*a) . 

Then the observable (S) can be written as 

1 1 dW(X,a 



A=1 

i 1 \ a dW 



Using eqs. (3.10) and (3.11), we obtain the one-loop result 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



The effective action can be obtained by restricting the diagrams that appear in the 
perturbative expansion of the free energy to one-particle irreducible (1PI) diagrams. It 
is the same as the free energy at the one-loop level, but starts to deviate at the two- 
loop level. In ref. [27] the two-loop effective action was calculated by evaluating the 1PI 
diagrams (a)~(d) in figure 1, and the two- loop effect turned out to vanish in the large- N 
limit. However, in order to calculate the observable (S), we need to evaluate the one- 
particle reducible (1PR) diagrams (f)~(h) in figure 1, which actually turn out to survive 
the large- N limit. 
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If we parametrize the higher order contributions Wj (j > 2) to the free energy in eq. 
(3.9) as 



Wj = -N 2 Wj(N) a^-fi 



the perturbative expansion of the observable (S) can be written as 
1 1 00 



(3.17) 



(3.18) 



J'=2 



The contribution of the 1PI diagrams to W2{N) is calculated as [27] 

w£ P1) (N) = ±{f 1 (N) + 4F 3 (N)} , 



(3.19) 



where the functions Fi(N) and F 3 (N) are defined by eqs. (A. 20) and (A. 22). The large- N 
behavior is found to be 



w 



(3.20) 



which vanishes in the large- N limit. 

On the other hand, the contribution 
of the 1PR diagrams to W2(N) are calcu- 
lated as 



w. 



(1PR), 



1 



2 (N) = ^(N) 



(3.21) 



cc 

Q. 



where F^N) is given explicitly by eq. 
(A.28). In figure 2 we plot ^ 1PR) (AT) 
against We find that 
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Figure 2: The coefficient u4 1PR ' ) (A^), which comes 
from the 1PR two-loop diagrams, is plotted against 
■jpz for N — 2, 3, 4, • • • , 32. The straight line repre- 



within the machine precision for TV = 
2, 3, 4, • • • , 32. Thus the 1PR diagrams 
yield the contribution to W2(N) which 
survives the large- N limit. Summing the 
contributions from the 1PI diagrams and the 1PR diagrams, we obtain the large- N behavior 



scnts w 



(1PR) 



(AO = 1 _ i 



W2 (N) = w { 2 m (N) + w { 2 1FR \N) = 1 + O ^ 



N 2 



(3.23) 



4. All order result from one-loop calculation 



In this section we apply a method in ref. [31] to obtain an all order result for the observable 
from one-loop calculation. The crucial point is that the free energy and the effective action 
are related to each other by the Legendre transformation. Therefore, one can obtain the 
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free energy by evaluating the effective action at its extremum. Since the effective action 
enjoys the one-loop dominance in the case at hand, we may obtain the free energy, and 
hence the observable, to all orders in J* in the large-TV limit. 

By expanding the theory around a rescaled fuzzy-sphere configuration = (3 L^, we 
obtain the one-loop effective action T in the large- N limit as 4 

1 ■ r lp-lap) + log/3, (4.1) 



lim -^T(P) 

n—>oc j\ z 



6 



where (3 = (3yN. The function of (3 on the right-hand side has a local minimum for 
a > from which one obtains the critical point d cr in eq. (2.8). The value of (3 which 

gives the local minimum can be obtained by solving a fourth order algebraic equation, and 
it can be written explicitly as 



= /(d) d = - d 



1 + VT+6+J2-6 + 



(4.2) 



where 



4d 3 



1 + A 1- 



512 

27 d 4 



+ 1-1/1- 



512 



27 a 



(4.3) 



Plugging this solution into the one-loop effective action (4.1), we obtain an all order result 
for the free energy as 



lim —pr W 

N->oo N 2 



/(« 



z\4 



1 



af(&r + log /(d) 



(4.4) 



(4.5) 



By using (3.15), we can readily obtain an all order result for the observable (S) as 

^< s > = ! - ii 5/ < 5 > 3 • 

where we have used the fact that (3 = f(a) extremizes the one-loop effective action (4.1). 

In order to check that the two-loop contribution obtained in section 3 can be reproduced 
correctly, let us expand the all order results at large a. First the solution (4.2) can be 
expanded in terms of 4j as 



/(«) 



a 



1 



oo 

i=i 



12 120 1456 



a 12 



a 



16 



The expansions for the free energy and the observable are obtained respectively as 

364 



1 jjr ^ , _ 1 14 HO 

lim — ; r W = — — + log a — — 7 — —77 — 



N^oo N 2 
1 



24 



a 



d 4 3d 8 3d 



12 



a 



16 



1 



28 



110 1456 



lim -={S) = -7T7 + 1 + ^J + ^ + ^W + ^ 



N^ooN 2 ' ' 24 d 4 3d 8 d 12 d 16 



+ • 



(4.6) 

(4.7) 
(4.8) 



The third term, which corresponds to the two-loop contribution, indeed agrees with the 
result (3.23), which we obtained by the direct calculation. 



'Here and henceforth we neglect an irrelevant constant term in the effective action and the free energy. 
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Let us discuss the convergence ra- 
dius of these series expansions. In figure 
3 we plot the ratio of the coefficients 
in (4.6), which is found to converge as 



lim 



9+1 



512 
~27 



(4.9) 
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This implies that the infinite series (4.6) 
converges for a > a cr . The expansion 
(4.7) for the free energy as well as that for 
the observable (4.8) has the same prop- 
erty. Thus in contrast to the perturba- 
tion theory in field theories, which usu- 
ally yields merely an asymptotic expan- 
sion, the perturbative expansion around the fuzzy sphere in the matrix model has a finite 
radius of convergence. The lower critical point (2.8) lies precisely on the convergence circle. 



Figure 3: The ratio of the coefficients in (4.6). 
The horizontal line represents ^ = 18.96296.... 



5. All order calculation of other observables 

So far we have focused on a particular observable (S), which can be obtained by differen- 
tiating the free energy. However, the method for deriving all order results is applicable to 
general observables [31] since the calculation reduces to the evaluation of the free energy 
with an appropriate source term. 

Suppose we want to calculate (O). Then we consider an action 

S e = S + eO, (5.1) 

and calculate the corresponding free energy W e , which has the expansion 

W e = W + e{Q) + 0(e 2 ) . (5.2) 

In order to calculate the free energy to all orders, we first calculate the effective action 
for the system (5.1) around the rescaled configuration = (3L„. Then the free energy 
W € can be obtained by evaluating the effective action at its extremum. Since the effective 
action enjoys the one-loop dominance, we obtain the free energy to all orders from the 
one-loop effective action. Let us expand the one-loop effective action F e as 

r e (/3) = r(^) + er 1 (^) + o(e 2 ) . (5.3) 

The "quantum solution" is given by solving the equation 

-~T e (P) = 0, (5.4) 
op 

whose solution is denoted as 

p = f(a) + eg{a) + 0{e 2 ) . (5.5) 
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Note that the first term is given by (4.2), and the second term represents a shift due to 
the source term. 

By plugging this solution into (5.3) and extracting the 0(e) term, we obtain the 0(e) 
term of the free energy, which is nothing but (O). Thus we obtain 



(5.6) 

/9=/(5) 



Note that the second term vanishes since (3 = f{a) extremizes T(/3), and the first term can 
be obtained by omitting the 1PR diagrams in the one-loop calculation of the observable 
(O) and replacing a by f(a). 

As a concrete observable let us consider the spacetime extent (-^tr(yl At ) 2 ). The one- 
loop result is given as [52] 

lim 1 / l t r(A„) 2 \ = - a 2 - \ . (5.7) 

The first (second) term corresponds to the classical (one-loop) contribution. As can be 
seen from (C.8) of ref. [52], the one-loop contribution is given totally by a tadpole diagram, 
which is one-particle reducible. Therefore, the all order result can be obtained as 

1/1 , . , 2 \ 1 ,,_, 2 1 ~2 1 5 48 572 

lim — ( — tr (A u 2 > = - f(a 2 = - a 2 5 = ^ tj . 5.8) 

N^ooN\N 1 4 Jl ' 4 a 2 a 6 a 10 a 14 y J 

Next we calculate the Chern-Simons term 

m = Yn w 1 < a » a » a p) ■ ( 5 - 9 ) 

The one-loop result in the large- N limit is given as [52] 

lim -^={M) 1 _ loop = -\ a 3 + 4 , (5.10) 



where the first (second) term corresponds to the classical (one-loop) contribution. Similarly 
to the case of (-^tr(A^) 2 ), the one-loop term comes solely from 1PR diagrams. Thus we 
obtain the all order result as 

1 , , x 1 „ 1 » 1 4 112 440 , 

lim —(M } = ~-f (a f = --a 3 + ^ + — + — ^ + ^3 +••• . 5.11 

The all order result for (■^ti(F IM/ ) 2 ) , where = i [A^A,,], can be readily obtained 
by using the exact result [52] 

ltr(F^) 2 ^+3a(M) = 3(l-^ . (5.12) 

Using (5.11), we obtain 
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From the relation 



w^ s) = \(^ F ^ 2 ) + 7N a{Mh (5 - 14) 

we find that (5.11) and (5.13) are consistent with the result (4.5) obtained in the previous 
section. Clearly the series expansions that appear in this section have the same radius of 
convergence as that for /(a). 



6. Comparison with Monte Carlo results 
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Figure 4: Monte Carlo results for various quantities are plotted against a for TV = 8, 16,32. The 
dotted (dashed) line represents the classical (one-loop) result, while the solid line represents the 
"all order" result. For the observable (jjtr(F^) 2 ), the one-loop result coincides with the classical 
result since there is no one- loop term; See eq. (5.13). 

In this section we compare the all order results (4.5), (5.8), (5.11) and (5.13) obtained 
in the previous section with our Monte Carlo data in ref. [52]. Figure 4 shows the Monte 
Carlo results for the four observables as a function of a for N = 8, 16, 32, where we also 
plot the classical, onedoop, and all order results. Note that the lines representing the all 
order results terminate at the critical point a = a cr . The all order results nicely reproduce 
the behavior near the critical point. This reinforces the validity of the all order calculation. 
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7. The dynamical gauge group revisited 



The 3d YMCS model (2.1) has var- 
ious classical solutions, which represent 
multi fuzzy spheres. Among them, the 
solutions 



A, 



(7.1) 




with N = nk, which correspond to k- 
coincident fuzzy spheres, are of particu- 
lar interest since they give rise to a non- 
commutative gauge theory on the fuzzy 
sphere with the gauge group of rank k 
[21]. If the true vacuum turns out to be 
described by the solution with some k, 
we may conclude that the gauge group of 
rank k has been generated dynamically. In ref. [52] we discussed this issue by comparing 
the one-loop effective action evaluated at the classical solutions (7.1). We found that the 
single fuzzy sphere always has the lowest free energy at the one-loop level in the fuzzy 
sphere phase. This implies that the dynamically generated gauge group is of rank one. 
From the view point of the all order calculation discussed in section 4, we should consider 
the rescaled configuration 



Figure 5: The all order result of the free energy 
for k coincident fuzzy spheres (k = 1,2, •••,6) is 
plotted against a. 



A li = /3LW®l k 



(7.2) 



and evaluate the one-loop effective action at its extremum, which gives the all order result 
for the free energy. In what follows we examine whether the higher order corrections to 
the free energy alter the conclusion. 

The one-loop effective action for the rescaled configuration (7.2) is given as [52] 



(7.3) 



Similarly to the single fuzzy sphere case (k = 1), the effective action (7.3) has a local 
minimum at 



P 



7 a ( 1 + y/l + + J2- <5« + - 2 

4 V V Vi+W) 



(7.4) 



where 



_4 

4d s 



1 + A/l 



512 k 2 



27a 



^4 



+ 1 



512 A; 2 
27 a 4 



(7.5) 



The local minimum exists if and only if a > 5$ = y 
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In figure 5 we plot the all order result for the free energy for k = 1,2, •••,6, which is 
obtained by evaluating the effective action (7.3) at its extremum (7.4). We find that the 
single fuzzy sphere (A; = 1) always has lower free energy than the coinciding fuzzy spheres 
(k > 2). Thus we conclude that the dynamically generated gauge group is of rank one even 
if we take account of the higher order contributions. This conclusion is consistent with our 
observations in the Monte Carlo simulation [52]. 



8. Summary 

In this paper we have clarified some peculiar aspects of the perturbative calculation around 
fuzzy sphere solutions in matrix models at large N. By direct two- loop calculation we have 
shown that the one-loop dominance, which holds for the effective action, does not hold 
for observables in general. However, we can obtain all order results for observables from 
onedoop calculation by shifting the center of expansion to the quantum solution, which 
extremizes the effective action. The validity of this method has been demonstrated by 
comparison with the direct twodoop calculation for the free energy and (S). We have also 
shown that the all order results for various observables reproduce nicely the behavior of 
our Monte Carlo data near the critical point. 

From the all order results we have found that the perturbative expansion around the 
fuzzy sphere has a finite radius of convergence. We recall that similar phenomena occur 
in exactly solvable matrix models in the planar large- N limit due to the exponential - 
rather than factorial — growth of the number of planar diagrams. In the present case the 
convergence property can be attributed to the onedoop dominance of the effective action. 

The conclusions listed above should also hold for other matrix models, which incorpo- 
rate higher-dimensional fuzzy manifolds. 
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A. Details of the two-loop calculation 

In this section we present the details of the two-loop calculation. In particular we explain 
how eqs. (3.19) and (3.21) are obtained. 

Let us note first that the kinetic term (3.6) can be written as 



S kin = iVtr 



(A.l) 
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where is an operator acting on a N x N matrix M as 

P /X M=[X IX ,M] (A.2) 

with being the fuzzy sphere solution (2.3). The operator (V\) 2 can be diagonalized by 
the so-called "matrix spherical harmonics" Y[ m (0 < I < N — 1, —l<m<l), which form 
a complete basis in the space of N x N matrices satisfying 

tr ( Y L Y l'™') = S W $mm' , (A.3) 

Yl = (-l) m Y^ m . (A.4) 
Similarly to the usual spherical harmonics, they also possess the properties such as 



m+l 



[L+,Y lm ] = v /(Z_ m )(/ + m + i)y Ii 

[£-,ljm] = V(J + m)(Z-m+l)yj, m _i , 

[L 3 ,Y lm ] =mY lm , (A.5) 

where L± = L\ ± zL2- From these properties we find that Y^ m is an eigenvector of the 
operator (V\) 2 , namely 

(Vx) 2 Y lm = a 2 l(l + l)Y lm . (A.6) 

By expanding the matrices , c and c as 

JV-l i N-l i N-l i 

Afj, = ^""^ A^i m Yi m , C = c lwXlm i C = ClmXlm j (A-7) 

Z=l m=—l Z=l m=—l 1=1 m=—l 

the propagators can be brought into the diagonal form 

1 (-l) m 

(^film^i/i'm')o = ]y^2 / + 1) VA^m-m' > ( A - 8 ) 
1 (-l) m 

(QmQ'm')o = jy^2 / -f 1) hj'^™-™' ' ( A - 9 ) 

where the symbol ( • )q represents a VEV using only the kinetic term Skm i n eq.(3.5). 
In the two-loop calculation, we use the identity 

= (-l) JV -V(2Zi + l)(2Z 2 + l)(2Z 3 + l) ( k h k )l ) ) ) ) , (A-10) 

\mi i7i2 rri3 I ILL 

where L is defined by JV = 2L + 1, and Wigner's (3j) and {6j} symbols are given explicitly 
[57] by Racah's formula as 

\ mi ni2 mz I 
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x - mi)! (h + mi)! (l 2 - m 2 )! (h + m 2 )\ (h ~ m 3 )l (Z 3 + m 3 )! 
x £(-!)*{*! (t - Z 2 + h + mi)! (* - Zi - m 2 + Z 3 )! (h +h~h- t)\ 



t 



h h h 
m\ m 2 m 3 



x(Zi-mi-t)!(Z 2 + m2-t)!} (A.ll) 

a/ A(lil 2 k) A(l 1 m 2 m 3 ) A(m 1 l 2 m 3 ) A(mim 2 l 3 ) 

£(-l)'(t + 1)! {(* - Zi - h ~ h)\ (t-h-m 2 - m 3 )\ 
t 

(t — mi — l 2 - m 3 )\ (t — mi — m 2 - Z 3 )i} 
| (h + h + mi + m 2 - t)\ (l 2 + l 3 + m 2 + m 3 — t)\ 
x(l 3 + h +m 3 + mi -t)l} \ (A.12) 



x 

X ' 
X 



where 



A(a6c) = (a + fe - c) / (b + C - a)! ^ + a - 6)! . (A.13) 
(a + b + c + ly. 

The sum of t is taken over all positive integers such that no factorial has a negative 
argument. If we have a negative argument in the factorial elsewhere, the (3j) and {6j} 
symbols are defined to be zero. From (A. 10) we also obtain the formula 

N-i h 

Ki,mii^2,m 2 ] = ^ ^2 fh,ml,l2,m^h,m3 , (A. 14) 

i 3 =l m3 =—l 3 

teh-m, = (-l)-3(_i)iv-i {1 _ (_i)'i+b-H.} > /(2I 1 + 1)(2Z 2 + 1)(2Z 3 + 1) 

(A.15) 



'/l,mi,/ 2 ,m2 

x (h h h \{h h h 

\ mi m 2 —m 3 ) | L L L 



With the help of these formulae, the contribution of the 1PI two- loop diagrams (a)~(d) 
to the coefficient w 2 (N) can be calculated as [27] 

4V) = -^Fi(N) , (A.16) 

w { 2 \N) = ^{Fi(N)-F 2 (N)} , (A.17) 

w®{N) = ^F 3 (N) , (A.18) 

4 d) (A0 = iFi(iV) , (A.19) 
where the functions Fi(N), F 2 (N) and F 3 (N) are given explicitly as 

JV-l 



FUN) = 1 V (2*i + l)(2*2 + l) 
U j iV,t- Zi(Zi + l)Z 2 (Z 2 + l) 

(1,12 = 1 
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U^i™,*^ '1^1 + + 1)^3 + 1) I J 

-\/ (Zi - mi)(Zi + mi + l)(l 2 + m 2 ){l2 - m 2 + 1) 

x / *i Z 2 /s V Zi h h \ 
Imi + 1 m 2 m 3 i I mi m 2 + 1 m 3 J 

+\V(h + mi)(h - mi + 1)(Z 2 - m 2 )(Z 2 + m 2 + 1) 

x / /i *2 is V ii ^2 h \ 
\ m\ — 1 m 2 m 3 / \ mi m 2 — 1 m 3 / 



+ mim 2 



' «i h h y 

mi m 2 m 3 / 



(A.21) 



F 3 (iV)- V fl ( D'i+fa-H.1 (2Zi + l)(2f 2 + 1)(2Z 3 + 1) \ h h hY . . 

Summing up the contributions (A.16)~(A.19) and using the identity Fi(N) = — 2F 2 (iV), 
we obtain eq. (3.19). 

Now let us calculate the contributions from the 1PR diagrams (f)~(h). We first eval- 
uate the tadpole obtained by contracting two A's in the second term in (3.7). It is given 
as 



T A = -iVajtr ([A^A v \[L^A v \)+tT (L4 M , A V }[L^A V }) +tr ([A^A V \[L^ A v \) 

N-l 

= -Not ^ ^([ Y hmi, Y l 2 m2}[ L ^^ Y hm 3 ]) 
h,h,h=l rni,m2,m 3 



^^,-^-lxl\m,\-^-vlim2^ylimz A v i imi Api 2m2 A u i., >mz + A^i imi A v i 2m2 A v i^ rn . i ) 



N-l 



a,f-, ^ ^ imi Z 2 (Z 2 + l) 

l ll l 2 =lm 1 ,m,2 y ' 

x{2 tr ([Y hmi , y zm2 ] [L„, E i2 _ m2 ]) + tr ([Y hm2 , Y h _ m2 ][L^ Y hmi })} . (A.23) 

The second term of (A.23) vanishes since 

= tr ([Lfj,, Yi imi [Yi 2Tn2 , Yj 2 _ m2 ]]) 
= tr {[L^Y hmi ][Yi 2m2 ,Yi 2 ^ m2 ]) 

+tr (Yi 1Tni [[L^i, Yi 2Jn2 \, Y\ 2 _ m2 \) + tr (Yi imi \Yi 2Jn2 , [L^, Yi 2 _ m2 ]]) 
= tr {[L^, Yi imi ]]Yi 2m2 , Yi 2 — m2 \) 

+tr(yi im j[L M ,y, 2tJW ],y, 2 _ m2 ])-tr(y, imi [[L M ,yi 2 _ m2 ],y, am2 ]) , (a.24) 



(_l)m 2 
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where the second and third terms in the last line cancel each other. Similarly we evaluate 
the tadpole obtained by contracting c and c in the fourth term in (3.7). It is given as 



N-l 



a N ^ ^ ^ ^ Q2m.2 j 4/i/im.i c /3m,3 tr {[L^, Y[ 2m2 \ [Y; imi , y^ 3m3 ]) 



2l,i2,23=l mi,m2,m 3 
N-l 



\rr12 



1 " ~ ( — l) r 

— y, y, Ani imi —— — — tr ([Lfj,, Yi 2 — m2 ]\Yi imi , Yi 2m2 ]\) 

h,l 2 =l m\,m,2 v ' 



(A.25) 



Thus we find that 



T gh = ~2 Ta ■ 



(A.26) 



By contracting two ^4's in the product of two tadpoles, we obtain the contribution 
from the diagrams (f)~(h) to the coefficient w 2 (N), which we denote as w^iN), w^(N) 
and w^\n), respectively. We obtain, for instance, 

N-l 



4\N) = 2o? £ £ 

h,te,h,l±=l mi,m2,m3,iri4 



(-l)™2+m 3 



Z 2 (*2 + l)fe(*3 + 1) 

tr ([^imu^2m 2 ][^! ^2-m 2 ]) tr ([^ 4 m 4 ) ^3m 3 ][^!/! ^-n^D 
2 ^_^jmi+rrt2+m,3 

JV,,^ , S Zi(Zi + l)Z 2 (/ 2 + l)/3(/3 + l) 

tr ([^liroi ) ^2m2][^i ^2— n^D tr (P'h — mi > ^ms] [-^/uj ^3—^13]) 

|^4(iV) , 



where we have defined 
F A (N) = 



N-l 

E E 

Zi,Z2,Z3=l "ii,m 2 ,m3 



(2Zi + 1)(2Z 2 + 1)(2Z 3 + l)(-l) m i+ m 2+ m 3 

Zl(/l + l)/ 2 (/ 2 + l)/ 3 (/3 + l) 



m 2 m 3 



^1 ^2 ^2 

mi m 2 — m 2 



h h h I 

L L L j 

h h h 

— 77l\ 7713 — m 3 



(h + m 2 )(^ 2 - m 2 + ~ m 3)(h + m 3 + 1) 



<*i ^2 ^2 
mi m 2 -m 2 + 1 



^1 ^3 ^3 
-mi m 3 — m 3 — 1 



+ 2" V - m 2 ){h +7712 + 1)(^3 + ™ 3 )(7 3 - m 3 + 1) 



il ^2 ^2 

mi m 2 — m 2 — 1 



^1 ^3 ^3 
-mi m 3 — m 3 + 1 



(A.27) 



(A.28) 
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Due to the relation (A. 26), the contributions from the other diagrams are readily obtained 
as 

4 g \N) = ^F 4 (N) , (A.29) 
= ~F*(N) ■ (A.30) 
Summing up all the contributions, we obtain eq. (3.21). 
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